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ABSTRACT 

The existence of ultra-fast rotating neutron stars (spin period P < 1 ms) is expected on the basis of current 
models for the secular evolution of interacting binaries, though they have not been detected yet. Their formation 
depends on the quantity of matter accreted by the neutron star which, in turn, is limited by the mechanism of mass 
ejection from the binary. An efficient mass ejection can avoid the formation of ultra-fast pulsars or their accretion 
induced collapse to a black hole. We propose that significant reductions of the mass-transfer rate may cause the 
switch-on of a radio pulsar phase, whose radiation pressure may be capable of ejecting out of the system most of 
the matter transferred by the companion. This can prevent, for long orbital periods and if a sufficiently fast spin 
has been reached, any further accretion, even if the original transfer rate is restored, thus limiting the minimum 
spin period attainable by the neutron star. We show that close systems (orbital periods P or b ~ 1 hr) are the only 
possible hosts for ultra-fast spinning neutron stars. This could explain why ultra-fast radio pulsars have not been 
detected so far, as the detection of pulsars with very short spin periods in close systems is hampered, in current 
radio surveys, by strong Doppler modulation and computational limitations. 

Subject headings: accretion, accretion discs - stars: binaries — stars: neutron — (stars:) pulsars: general — 
X-ray: stars — X-ray: general 



1. INTRODUCTION 

The shortest spin period ever observed for a rotating neu- 
tron star (hereafter NS), P min = 1.56 ms (Backer et al. 1982), 
is not a physical limit for the stability of spun-up NSs. In fact, 
the period P m ; n is longer than the limiting period, Pij m , below 
which the star becomes unstable to mass shedding at its equator. 
Cook, Shapiro & Teukolsky (1994) have shown that f\ im — 0.5 
ms for most equations of state (EoS) adopted to describe the 
ultra-dense nuclear matter. 

The re-acceleration of a NS to ultra-short periods depends 
remarkably on the amount of mass (and hence of angular mo- 
mentum) accreted. Typically M\ ms ~ 0.35 M Q must be accreted 
to reach P ~ 1 ms (e.g. Burderi et al. 1999). Most donor stars 
in systems hosting recycled millisecond radio pulsars (MSP) 
have certainly lost, during their interacting binary evolution, 
a mass ;> M\ ms . They now appear as white dwarfs of mass 
~ 0.15-0.30 M (e.g. Taam, King & Ritter 2000), whose pro- 
genitors are likely to have been stars of ~ 1.0-2.0 M Q (Web- 
bink, Rappaport & Savonije 1983; Burderi, King & Wynn 1996, 
Tauris & Savonije 1999). Examples of possible evolutionary se- 
quences, in which the assumption of conservative mass transfer 
is plausible, are computed in §2. Hence, in order to explain 
the lack of ultra-fast rotations, we have to find physical mecha- 
nisms able to prevent the accretion of a considerable fraction of 
the mass lost by the companion (M; osl = 0.8 — 1.7 M ) onto the 
NS. 

In this Letter we propose a mass ejection mechanism based 



on the sweeping effects of the energy outflow from a rapidly 
spinning NS undergoing a radio pulsar phase. The sweeping 
effects of the energy outflow from a rapidly spinning NS have 
been firmly established since the early works of Shvartsman 
(1970) and Illarionov & Sunyaev (1975). Ruderman, Shaham 
& Tavani (1989) and Shaham & Tavani (1991) discussed in par- 
ticular the case of recycled pulsars in low mass binary systems. 
As soon as the accreting plasma moves out beyond the light- 
cylinder radius (where an object corotating with the NS attains 
the speed of light, R^q = cP/2tt), the NS becomes generator of 
magnetodipole radiation and relativistic particles, whose pres- 
sure may expel the matter overflowing the Roche lobe (see also 
Campana et al. 1998). In the following we determine the de- 
pendence of this mass ejection mechanism on the parameters of 
the system (§3) and show that it naturally provides both an ex- 
planation for the values of the mass and rotation of the observed 
MSPs (Thorsett & Chakrabarty 1999; Tauris & Savonije 1999) 
and an indication about where ultra-fast spinning NSs could re- 
side (§4), suggesting in turn why they could have been elusive 
up to now. 

2. ACCRETION AND EJECTION DURING MASS TRANSFER 

It is well known that an upper limit to the accretion rate 
is given by the Eddington limit. However, for typical initial 
masses ,$,1.6 M Q and initial orbital periods P or b ;$ 10 days, the 
donor transfers mass at iM^-Eddington rates, in principle mak- 
ing all Mi os, available for the recycling of the NS. We explored 
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the conservative mass transfer scenario computing the system 
evolution (with initial and final parameters listed in Table [j]) 
with the ATON1.2 code (D'Antona, Mazzitelli, & Ritter 1989). 
The mass loss rate is computed following the formulation by 
Ritter (1988), as an exponential function of the distance of the 
stellar radius to the Roche lobe, in units of the pressure scale 
height. This method also allows to compute the first phases 
of mass transfer, during which the rate reaches values which 
can be much larger than the stationary values, due to the ther- 
mal response of the star to mass loss. For systems having a 
mass donor of 1 .2M Q we assumed either that the donor fills the 
Roche lobe while it is evolving towards the red giant branch 
(case B) or during the core-hydrogen burning phase (case A). 
Results are shown in Table 1 . 

In cases 2, 3 and 4, mass transfer is guided by the thermal 
and nuclear evolution of the secondary; in case 1 it is guided 
by magnetic braking. The latter evolution is expected to stop at 
f or b ~ 2.5 hr, when the star becomes fully convective. It may 
resume at a shorter period guided by gravitational wave loss 
of angular momentum, like in the cataclysmic binary systems, 
reaching a minimum period of ~ 1.03 hr. For the sequences 
2, 3 and 4, the final period is in the range 90 - 150 hr; the 
final mass of the donor is in the range 0.21 -0.33 M Q (see 
Tab. [[]), in agreement with the masses of the companion stars 
of MSPs as inferred from accurate timing measures (e.g. Bur- 
den, King, & Wynn 1998). On the other hand the assumption 
of conservative mass transfer leads to final NS masses in the 
range 2.3-2.7 M , much larger than the inferred masses of the 
NSs in these systems and dangerously close to (or even larger 
than) the maximum mass allowed for a NS in most of the pro- 
posed EoS (Cook, Shapiro & Teukolsky 1994). This suggests 
that either the mass transfer cannot be conservative or that the 
final NS must be very massive implying that accretion-induced 
collapse to a black hole is a likely outcome for LMXBs. Indeed 
the masses of a sample of radio pulsars in binary systems were 
estimated by Thorsett & Chakrabarty (1999). These are con- 
sistent with a remarkably narrow Gaussian distribution, with 
Mns = 1.35 ±0.04 Mq, although the sample is contaminated by 
five relativistic NS-NS binary systems, the progenitors of which 
are massive X-ray binaries. If the proposed value of the amount 
of mass accreted, < 0.1 Mq, is representative of the MSP ob- 
served sample, this suggests that significant mass losses occur 
during the Roche lobe overflow phase. 

A propeller effect (Illarionov & Sunyaev 1975) has been of- 
ten invoked to explain this discrepancy. In the widely accepted 
scenario of accretion onto a magnetized NS (see e.g. Hayakawa 
1985 for a review) the accretion disc is truncated at the magne- 
tospheric radius R m , at which the magnetic field pressure equals 
the pressure of the matter in the disc. The value of R m for a 
Shakura-Sunyaev accretion disc (e.g. Burderi et al. 1998) is, in 
a first approximation, close to ~ 0.5 times the Alfven radius R^: 

R m ~ 0.5 X R A = 1.1 X 10V26 7 ^6 2/7L 37 /7ml/7cm ' Where L 37 is 

the accretion luminosity in units of 10 37 erg/s, m is the NS mass 
in solar masses, R& is the NS radius in units of 10 6 cm, /j,26 is the 
magnetic moment of the NS in units of 10 26 G cm 3 (/j, = B S R 3 
with R and B s NS radius and surface magnetic field along the 
magnetic axis, respectively). Accretion onto a spinning mag- 
netized NS is centrifugally inhibited once R m lies outside the 
corotation radius Rqo, the radius at which the Keplerian an- 
gular frequency of the orbiting matter is equal to the NS spin: 
Rco = 1-50 x \0 6 m l l 3 P^ cm, where P_3 is the spin period in 
milliseconds. In this case a significant fraction of the accret- 



ing matter could in principle be ejected from the system: this is 
called propeller phase. 

The virial theorem sets stringent limits on the fraction of mat- 
ter that can be ejected in this phase. In fact, it states that, at 
any radius in the disc, the virialized matter has already liber- 
ated (via electromagnetic radiation) half of its available energy. 
Considering that R m ~ /?ns for MSPs, the matter at the mag- 
netosphere has radiated ~ 50% of the whole specific energy, 
£acc = GM/Rns, obtainable from accretion. To eject the same 
matter (close to the NS surface) ^1/2 £ acc must be given back 
to it. As the only source of energy is that stored in the NS 
by the accretion process itself, the typical ejection efficiency is 
~ 50%. Actually, with a fine tuned alternation of accretion and 
propeller phases of the right duration, the ejection efficiency 
can be higher than 50%, although a difficulty with this scenario 
is that, once the system has reached the spin-equilibrium, no 
further spin up takes place and the storage of accretion energy 
in a form that allows its subsequent re-usage for ejection is 
impossible. Thus the accreted mass is about a half of the trans- 
ferred mass, i.e. ~ 0.4-0.8 Mq. This has two main conse- 
quences: a) even taking a propeller phase into account, NSs in 
MSPs will either be very massive or even collapse into black 
holes, b) as the amount of mass accreted is considerably larger 
than the minimum required to spin up the NSs to ultra short pe- 
riods, one has to invoke an ad hoc final, long lasting propeller 
phase with a highly effective spin down to form the observed 
population of moderately fast spinning MSPs. 

An alternative viable hypothesis to explain the lack of ultra- 
fast rotating NS is that gravitational wave emission balances the 
torque due to accretion (see a review in Ushomirsky, Bildsten, 
& Cutler 2000). However, while these emission mechanisms, 
whose importance in LMXBs will be probably clarified in the 
next future by detectors such as LIGO-II, might explain the ob- 
served spin periods of spun-up NSs, which seem to cluster be- 
tween 260 and 590 Hz (see e.g. van der Klis 2000), they cannot 
solve the problem of the large final masses of NSs. 

The only way to overcome these difficulties is to obtain ejec- 
tion efficiencies close to unity. This is indeed possible if the 
matter is ejected so far away from the NS surface that it has an 
almost negligible binding energy GM/r <C E- dcc . As the NS is 
spinning very fast, the switch on of a radio pulsar is unavoid- 
able once R m > R^c- The pressure exerted by the radiation field 
of the radio pulsar may overcome the pressure of the accre- 
tion disc, thus determining ejection of matter from the system. 
Once the disc has been swept away, the radiation pressure stops 
the infalling matter as it overflows the inner Lagrangian point, 
where £„ ~ 0. 

3. THE EFFECT OF THE PULSAR ENERGY OUTFLOW 

The push on the accretion flow exerted by the (assumed 
dipolar) magnetic field of the NS can be described in terms 
of an outward pressure (we use the expressions outward or 
inward pressures to indicate the direction of the force with 
respect to the radial direction): Pmag = B 2 /An = 7.96 x 
10 14 /U26 r 6 6 dy/cm 2 , where r^ is the distance from the NS cen- 
ter in units of 10 5 cm. If the disk terminates outside 7?lc, 
the outward pressure is the radiation pressure of the rotat- 
ing magnetic dipole, which, assuming isotropic emission, is 
Ppsr = 2.04 x 10 12 P^^ 6 rg 2 dy/cm 2 . In Figure 1 the two out- 
ward pressures (Pmag and Ppsr) are shown as bold lines for 
typical values of the parameters (see figure caption). 

The flow, in turn, exerts an inward pressure on 
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the field. For a Shakura-Sunyaev accretion disc (see 
e.g. Frank, King & Raine 1992): P DISK = 1.02 x 

lO^a^/^no.^s^T 720 ^ 1740 ^" 720 / 1775 ^ 2178 dy/cm 2 , where a 
is the Shakura-Sunyaev viscosity parameter, «o.6i5 = n/0.615 ~ 
1 for a gas with solar abundances (where n is the mean particle 
mass in units of the proton mass), and / = [1 -{Rb/ r^) 1 ' 1 ] 1 ^ < 
1. As before we measure M in units of L37, supposing that 
all the matter reaches the NS surface. Note that this equation 
is valid in the zone C of the dik, that is at radii larger than 
3.45 x 10 7 cm for a luminosity of 10 37 erg s" 1 (Burderi, King, 
& Szuszkiewicz 1998). In Figure 1 the inward disc pressure for 
a luminosity Lmax, corresponding to a high accretion rate, is 
shown as a thin solid line. The disc pressure line, which inter- 
sects Pmag at Rlc, defines a critical mass-transfer rate M sw i tC h 
(i.e. a critical luminosity L sw j tc h, shown as a dashed line in 
Fig. 1) at which the radio pulsar switches-on. 

The intersections of the Pdisc line corresponding to Lmax 
with each of the outward pressure lines define equilibrium 
points between the inward and outward pressures. The equi- 
librium is stable at r = R m , and unstable at r = Pstop, which can 
be derived equating P P sr and Pdisk: 



W ~ 8 x 10 1 V 36 / 25 ^* 4/25 / 136/25 
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cm. 
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In fact, as Pmag is steeper than Pdisc, if a small fluctuation 
forces the inner rim of the disc inward (outward), in a region 
where the magnetic pressure is greater (smaller) than the disc 
pressure, this results in a net force that pushes the disc back 
to its original location R m . As Ppsr is flatter than Pdisc, with 
the same argument is easy to see that no stable equilibrium is 
possible at /?stop and the disc is swept away by the radiation 
pressure. This means that, for r > Rstop, no disc can exist for 
any luminosity < Lmax- 

It is convenient to divide the systems into "compact" and 
"wide" depending whether the primary Roche-lobe radius (Rl\) 
lies inward or outward Rstop, as they behave very differently in 
response to significant variations of the mass-transfer rate. The 
dependence of Rn on the orbital parameters is given by: 



R 



LI 



= 3.5 x lO w F? /3 (m+m 2 ) l/3 



1-0.462 



m 2 



m + m.2 



1/3 



using the approximation given by Paczynski (1971), where P^ 
is the orbital period in hours and to, TO2 are the NS and the com- 
panion masses in M Q , respectively. Therefore, the P or b that sep- 
arates "compact" and "wide" systems can be obtained imposing 
^stop/^li = 1 and solving forP orb : 

P crit = 1.05 x (a^sRff l50 ^rn^^ l5 PT 



1-0.462 



m.2 

TO+TO2 



1/3 



-3/2 



(to + to 2 )" 1/2 h (2) 



where L36 is Lmax in units of 10 erg/s. (In the following, al- 
ways keep in mind that the separation between "compact" and 

"wide" systems depends on P^ 5 , and has no 'absolute' mean- 
ing). 

When the luminosity alternates between its maximum value 
Lmax and a minimum luminosity Lviin < L sw ; tc h the behavior 
of a "compact" system is cyclic. During the high state, the 



magnetospheric radius is smaller than both the corotation and 
the light-cylinder radius and the NS will normally accrete mat- 
ter and angular momentum, thus increasing its spin (accretion 
phase). The sudden drop in the mass-transfer rate to Lviin ini- 
tiates a phase, that we termed "radio ejection", in which the 
mechanism that drives mass overflow through L\ is still active, 
while the pulsar radiation pressure prevents mass accretion. As 
the matter released from the companion cannot accrete, it is 
now ejected as soon as it enters the Roche lobe of the primary. 
When the system goes back to the high state, the disk might 
form again and the accretion can resume. 

The response of a "wide" system to the same variations of the 
accretion rate is quite different. When the mass-transfer rate re- 
covers the value corresponding to the high state, the cyclic be- 
havior is lost. Indeed, since in this case R^\ is located beyond 
^stop, ^toisc < ^psr at even in the high state and the accre- 
tion cannot resume. This means that for a "wide" system, once 
a drop of the mass-transfer rate has started the radio ejection, a 
subsequent restoration of the original mass-transfer rate is un- 
able to quench the ejection process (as already pointed out by 
Ruderman, Shaham, & Tavani 1989). 

In conclusion, while the evolution without a radio ejection 
phase implies that a large fraction of the transferred mass is ac- 
creted onto the NS (because of the constraints imposed by the 
virial theorem), we have demonstrated that the switch-on of a 
radio pulsar (associated to a significant drop in mass transfer) 
could determine ejection efficiencies close to 100%, as the mat- 
ter is ejected before it falls into the deep gravitational potential 
well of the primary. For "compact" systems, we have shown 
that radio ejection is swiftly quenched by a resumption of the 
original mass-transfer rate, leading to the prediction that the 
amount of mass accreted is substantial. On the other hand, if 
a radio ejection starts in a "wide" system, this implies that the 
accretion is inhibited in the subsequent evolution. 

4. WHERE TO SEARCH ULTRA-FAST SPINNING NS 

A statistical analysis based on the current samples of detected 
MSPs (Cordes & Chernoff 1997) proved that, using different 
hypotheses for the period distributions of these sources, there is 
always a non negligible probability for periods P < P m i„. Pos- 
senti et al. (1999, 1998) performed population synthesis cal- 
culations including propeller and randomly choosing the mass 
accreted onto the NS in the interval 0.01 -0.4 M Q (this corre- 
sponds to efficiency of ejection between 60% and 99%). They 
confirmed that the process of recycling in low-mass binaries 
can produce a significant amount of ultra-rapidly rotating ob- 
jects, under different assumptions for the evolution of the mass 
transfer and the magnetic field. However, despite the large ef- 
forts devoted in the last years (D' Amico 2000; Crawford, Kaspi 
& Bell 2000; Edwards et al. 2001) no pulsar with P < P min has 
been observed so far. 

The more rapid the NS spins, the smaller is the drop of M 
needed to switch the pulsar on and expel the accretion disc (see 
§3). Once this occurred a new phase of accretion is possible 
only if the orbital period is short enough. Eq. (^|) predicts that 
spinning a NS up to P 1 ms requires very close (P or b ^1.5 
hrs) binary systems. 

The orbital Doppler shift on the pulsar signal can provide a 
natural observational bias against the detection of ultra-rapidly 
spinning pulsars. All the codes for searching pulsations from 
a source in a close binary system are a compromise between 
computational capability and sensitivity: on each data set, they 
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must perform a two dimensional search in the space of the un- 
known parameters dispersion measure (DM, related to the dis- 
tance of the object) and acceleration (resulting from binary mo- 
tion). Short data sets reduce the Doppler period shifts during 
the observation and relax the computational requests, but at the 
price of limiting the sensitivity (see e.g. Camilo et al. 2000). As 
a consequence, up to now orbital periods shorter than ~ 90 min 
have been poorly searched, even in the more favorable case of 
targeted searches (as those pointing to globular clusters), where 
one of the two parameters (DM) is known. This observational 
bias could be enhanced in presence of eclipses (favoured in very 
close binary pulsar systems, Nice et al. 2000) or in case of 
a large duty cycle of the pulsar signal, a low radio luminosity 
and a strong interstellar scintillation, already suggested for the 
elusiveness of ultra fast rotating pulsars (Possenti 2000). 

Even in the favorable case P or t, < P cl ; t , a steady accretion dur- 
ing all the Roche lobe overflow phase (i.e. no ejection episodes) 
could prevent the formation of a ultra-fast spinning NSs: in fact 
if the EoS is soft the NS could undergo a collapse to black hole. 



Actually, Cook, Shapiro & Teukolsky (1994) showed that (be- 
cause of the centrifugal pull determined by the rapid NS rota- 
tion) a NS can accrete matter even beyond the limit at which 
the black-hole collapse takes place. However, when the accre- 
tion halts and the fast spinning NS loses rotational energy via 
magnetodipole emission, a radial instability can set in leading 
to the formation of a black-hole (see Stella & Vietri 1999). 

In summary we have shown that: 1 . For very fast rotating 
NSs, if variabilities in the matter flow trigger an extended radio- 
pulsar phase, accretion onto a NS cannot recover unless also the 
orbital period is quite short. Thus extremely compact binary 
systems are strongly favoured for harbouring ultra fast rotating 
pulsars. 2. If the equation of state for the nuclear matter is soft 
enough, binary evolution with steady accretion rate would often 
lead to black-hole formation, via accretion induced collapse. 

This work was partially supported by a grant from the Italian 
Ministry of University and Research (Cofin-99-02-02). 
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Table 1 

Mass Evolution in Low Mass Binaries. 



System type M NS (M Q ) M donor (M Q ) P or b(hi&) 

Initial 1 1.40 1.20 10.46 

Final 2.58 0.02 1.42 

Initial 2 1.40 1.20 30.55 

Final 2.35 0.25 338.1 

Initial 3 1.40 1.40 33.95 

Final 2.54 0.26 483.8 

Initial 4 1.40 1.60 34.00 

Final 2.67 0.33 346.0 
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